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We describe a proposed experimental search for exotic spin-coupled interactions using a solid- 
state paramagnetic insulator. The experiment is sensitive to the net magnetization induced by the 
exotic interaction between the unpaired insulator electrons with a dense, non-magnetic mass in close 
proximity. An existing experiment has been used to set limits on the electric dipole moment of the 
electron by probing the magnetization induced in a cryogenic gadolinium gallium garnet sample on 
application of a strong electric field. With suitable additions, including a movable source mass, this 
experiment can be used to explore “monopole-dipole” forces on polarized electrons with unique or 
unprecedented sensitivity. The solid-state, non-magnetic construction, combined with the low-noise 
conditions and extremely sensitive magnetometry available at cryogenic temperatures could lead to 
a sensitivity over ten orders of magnitude greater than exiting limits in the range below 1 mm. 

PACS numbers: 32.10.Dk, ll.30.Er, 77.22.-d, 14.80.Va 


Introduction. — Experimental searches for macroscopic 
forces beyond gravity and electromagnetism have re¬ 
ceived a great deal of attention in the past two decades. 
Present limits allow for unobserved forces several mil¬ 
lion times stronger than gravity acting over distances of 
a few microns. Predictions of unobserved forces in this 
range have arisen in several contexts, including attempts 
to describe gravity and the other fundamental interac¬ 
tions in the same theoretical framework. For compre¬ 
hensive reviews, see [H-Q. The sub-millimeter range has 
been the subject of active theoretical investigation, no¬ 
tably on account of the prediction of “large” extra di¬ 
mensions at this scale which could explain the hierarchy 
problem Q. The fact that the dark energy density, of 
order (1 meV)^, corresponds to a length scale of about 
100 ^m also encourages searches for unobserved phenom¬ 
ena at this scale. Many theories beyond the Standard 
Model possess extended symmetries that, when broken 
at high energy scales, give rise to light bosons with ve^ 
weak couplings to matter. Examples include moduli 
dilatons [6| , and the axion-a light pseudoscalar motivated 
by the strong CP problem of QCD Q. These particles 
can generate weak, relatively long-range interactions be¬ 
tween samples of ordinary matter, including interactions 
that couple to spin. In a seminal paper Q, Moody and 
Wilczek derived three possible interactions for the ax- 
ion, and proposed searches sensitive to the T-violating 
“monopole-dipole” interaction between polarized and un¬ 
polarized test masses. 

We propose an experimental search for exotic spin- 
coupled interactions using a solid-state paramagnetic in¬ 
sulator as a detector. The candidate material is gadolin- 
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ium gallium garnet (Gd 3 Ga 50 i 2 , or GGG), which has 
been used as a detector in an experimental search for 
the electric dipole moment (EDM) of the electron Q. 
(See [i3,[nl for realizations of other EDM materials.) In 
that experiment, the signal is an induced sample magne¬ 
tization on application of a strong external electric field; 
in our proposal the magnetization is induced by an exotic 
monopole-dipole interaction with a dense, non-magnetic 
mass brought into close proximity. An exotic field cou¬ 
pling to the electron spins of an atom or an ion with spin 
S = \/s(s + 1) in the solid sample leads to a net spin 
excess in the sample. At the location of a particular ion 
in the sample, the spin excess ratio is given by [Hill 
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where mg is the magnetic quantum number, u is the lo¬ 
cal exotic field energy per ion, T is the sample temper¬ 
ature and k is Boltzmann’s constant. In our proposal, 
the energy shift predicted from the exotic coupling at 
the current experimental limits is much smaller than the 
thermal energy. However, the cumulative effect from the 
large number of electrons in a macroscopic, solid-state 
sample leads to a net spin alignment on the order of 10® 
Bohr magnetons cm“®. The resulting sample magneti¬ 
zation can be probed with great sensitivity using sensors 
based on superconducting quantum interference device 
(SQUID) technology. 

The idea of using solid-state materials with bound but 
unpaired electron spins to probe exotic fields was first 
proposed by Shapiro in 1968 [l^ in the context of EDM 
searches. Recently, other fundamental applications of 
these materials have been proposed, including tests of 
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Lorentz and CPT invariance |15l| . Experiments using 
EDM techniques, but on a different class of solid-state 
materials, have been proposed to search for cosmic ax- 
ions [l^ . Many experiments have been performed to 
search for spin-dependent macroscopic interactions us¬ 
ing other methods. Examples include NMR-type exper¬ 
iments sensitive to precession frequency shifts in various 
materials, including the paramagnetic salt TbFa [13, Hg 
and Cs comagnetometers [l^, polarized ^^®Xe and ^^^Xe 
gas [l3|, and polarized ^He gas [20l - l^ . in the presence 
of polarized and unpolarized masses. Other exp eriments 
search for effects in torsion pendulums |23l - l^ . neutron 
bound states in the Earth’s gravitational field [13, |13 , 
and longitudinal and transverse spin relaxation of polar¬ 
ized neutrons and ^He [l^-js^. An overview can be found 
in [ 3 . Other parameters being equal, the proposed solid- 
state technique affords an enhancement factor on the or¬ 
der of the Avogadro number relative to experiments in 
dilute vapor systems, though the latter are primarily sen¬ 
sitive to polarized nucleon couplings. With suitable con¬ 
trol of systematic effects, the ultimate sensitivity is more 
than 10 orders of magnitude greater than current labo¬ 
ratory limits in the range below 1 cm, and the technique 
is sensitive to exotic interactions of electrons presently 
unconstrained by either laboratory experiments or astro- 
physical observations. 


A study by Dobrescu and Mocioiu [3^ of the pos¬ 
sible interactions between non-relativistic fermions as¬ 
suming only rotational invariance revealed 15 forms 
for the potential involving the fermion spins. Nine of 
these are spin-spin interactions, which would necessi¬ 
tate spin-polarized test masses with low intrinsic mag¬ 
netism 1^ ; here we concentrate on monopole-dipole 
interactions between polarized and unpolarized objects. 
In the zero-momentum transfer limit, the possible inter¬ 
actions between a polarized electron and an unpolarized 
atom or molecule of atomic number Z and mass A are 
(in SI units, and adopting the numbering scheme in [13): 
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Here S = ‘ha/2 is the electron spin, h is Planck’s con¬ 
stant, f = r/r is a unit vector along the direction be¬ 
tween the electron and atom, v is their relative velocity, 
c is the speed of light in vacuum, me is the electron mass, 
and A is the interaction range. The factors gp and are 
the electron pseudoscalar and axial vector coupling con¬ 
stants, and gg and gy are the nucleon scalar and vector 
coup ling s. The couplings in Eq. [2]are not the most gen¬ 
eral we have included those for which the proposed 
experiment will likely have the greatest discovery poten¬ 
tial. We note that V 9+10 can also proceed via a spin-1 
interaction, in which case the coupling is g%gy ■ V4-1-5 can 



FIG. 1. Schematic of proposed experiment. A dense, non¬ 
magnetic planar source mass is modulated vertically above a 
paramagnetic detector mass, inducing a magnetization read 
out with the pickup coil and SQUID sensor. The detector is 
kept at temperatures near 1 K. 


also proceed via a spin -0 interaction, in which case the 
coupling is gggg (and the expression for V 4+5 above is 
scaled by A/Z). 

The experiment is illustrated schematically in Fig. [T] 
It is based directly on the apparatus used in Q; many 
parameters have been retained for the purposes of design¬ 
ing a practical device. A solid paramagnetic insulating 
sample or “detector” mass, in the form of a block or disk, 
is mounted in the sample space of a large dilution refrig¬ 
erator. A dense, unpolarized, non-magnetic insulator of 
similar size and shape is brought into close proximity and 
serves as “source” mass. The source-detector gap can 
be modulated (e.g., via translation or rotation stages) 
from 0-1 cm with the resolution of a few microns. The 
non-magnetic design eliminates the need for shielding be¬ 
tween the test masses. At closest approach, the source is 
essentially in contact with the detector, permitting force 
searches with potentially unprecedented sensitivity in the 
range below 1 mm. The detector magnetization induced 
by the exotic interaction with the source is sensed by a 
pickup coil surrounding the detector. The coil is coupled 
to a sensitive SQUID magnetometer. 

For the detector material we assume GGG, as in [3 ■ 
The use of GGG was first proposed in [33 and realized 
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in [^ . The material is chosen to maximize the induced 
spin density. GGG contains a high density of Gd^+ ions 
{fv 10^^/cm^), each of which has seven electrons in the 
4f shell: the shell is thus half-filled and all electrons are 
unpaired [1^. This property leads to a relatively strong 
magnetic response to external fields. For the source mass 
we assume bismuth germanate (Bi 4 Ge.'^Oi 2 , or BGO), a 
high density, non-magnetic insulator [22j| . 

Experimental Sensitivity .— An approximate sensitiv¬ 
ity and figure-of-merit can be derived using a simplified 
geometry; we examine the case of the V 9+10 interaction, 
which is the most widely studied. The interaction energy 
between an unpolarized, flat-plate source mass parallel 
to a flat-plate paramagnetic detector mass is given by: 
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where nm is the number density of molecules in the source 
mass and is their mass number, Ud is the number 
density of polarized ions in the detector mass and Sd is 
the ion spin, ad is the detector area, z{t) is the source- 
detector gap, and t^ and td are the source and detector 
thicknesses. This expression is exact for the case of a 
source of infinite area, in which case the detector spins 
align toward the source, in the direction normal to the de¬ 
tector plane. For the practical case of a finite source (con¬ 
sidered in detail below), there are edge corrections (which 
are the dominant effect in the 14-1-5 case; see Eq. ID). The 
proposed experiment is sensitive to changes in the in¬ 
duced field as the source mass is modulated, which we 
assume to occur at some fixed frequency oj significantly 
above the 1/f noise corner of the SQUID. The noise cor¬ 
ner was about 0.1 Hz for the SQUID in [^; we assume 
uj/2tt = 1 Hz. From Fourier analysis, the amplitude of 
the energy change per detector ion at lu is 

j-2 
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where wg-i-io = t/g+io/A^d and Nd = Udadtd is the num¬ 
ber of detector spins. In Eq. |4] we have taken the limit 
X ^ td, so that iig+io is a good approximation at the 
location of any particular ion in the detector mass (or, 
equivalently, in a layer of thickness ~ A at the top sur¬ 
face of a thicker detector). Here, /i(zo/A) is the modified 
Bessel function, and we have used z(t) = z + ^osinwt 
where z = ( 2 max + 2 :min )/2 is the average source-detector 
gap and zq = (zmax — is the amplitude of the 

source mass motion. 

In the limit ug+io kT, the spin excess ratio (Eq. [J) 
can be approximated as R = ug+io/(3fcT), whereby the 
total spin excess in the layer is SdNdUg+io/{ikT). The 
layer magnetization normal to the plane is 

99'BSdndU9+io ..X 
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where 0 = 2 for electrons and fip '^s the Bohr magneton. 
As in 0 , we assume the plane of the pickup coil to be sit¬ 
uated just below the bottom surface of the detector mass 
(Fig. [T|). The field of the magnetized slab of thickness A 
just below the bottom is slowly-varying across the sur¬ 
face, with a minimum at the center of ~ fioMX/^/od- 
Assuming a coil area of , a conservative estimate of the 
induced flux through the coil is thus 4)^ r; pgfMX^/od, 
where we include a suppression factor / for sub-optimal 
coupling between the detector and coil. 

In an experiment using a practical SQUID sensor, the 
pickup coil connects to a built-in input coil on the sensor. 
The changing flux from the detector induces a current, 
which flows into the input coil and produces a flux that 
couples inductively to the SQUID loop and is transduced 
to a voltage. The relationship between the magnetic flux 
picked up from the detector mass and the induced flux 
in the SQUID loop is @ 

M 
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where Lp and Li are the inductances of the pickup and 
input coils, respectively, and M is the mutual inductance 
between the input coil and SQUID. The factor (3 is the 
coupling efficiency which quantifies the flux reduction 
when 4>i is delivered to the SQUID sensor. 

The sensitivity of the experiment is based on the expec¬ 
tation that essentially all backgrounds can be suppressed 
below the intrinsic noise of the SQUID sensor. In an ex¬ 
periment limited by this noise (</>„, expressed in terms 
of magnetic flux per root bandwidth), the signal-to-noise 
ratio is given by: 


SNR=^^^, (7) 
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where t is the integration time. The sensitivity is calcu¬ 
lated by setting SNR = 1 and solving for gpgs- Gom- 
bining Eqs. [3] through [3 the result is: 
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Here we have used nm = Pm/ {oAm), where pm is the 
mass density of the source and u is the atomic mass unit. 
The factor Xd is the effective susceptibility of the detector 
mass: 
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and we have used Xd —t Ad/(1 + L>Xd) in Eq. |8l where 
D is the demagnetization factor. The efficiency e = 
/i(zo/A)e“^/'^(l — in Eq. [5] is of order unity for 

an optimized vertical geometry. 

A plot of g%gg vs A from Eq. [8] is shown in Fig. [3l 
using the parameters in Tables U and m We assume the 
source and detector can be brought into contact but also 
that these elements can only be made flat to within a 
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few microns, thus the minimum gap is set to 10 ^m. For 
each value of A, the maximum separation Zmax is chosen 
to maximize e, which occurs around Zmax = 3A. The 
efficiency is maximized for tm ^ oo, however, since very 
little is gained for tm > A, we set = 1 cm, about twice 
the maximum range of interest. A suppression factor of 
/ = 0.4 has been used, in accordance with Q. For a 
rectangular prism of square cross-section ad magnetized 
parallel to the thickness (here ~ A), we have used the 
approximation D Ri ^/^/+ 2A) [40l |. 

From Eq. |S1 an approximate figure-of-merit for the ex¬ 
periment is: 
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illustrating the importance of high pm, and z ~ zq ^ \ 
(through e). Over the range of interest, D varies be¬ 
tween 0.75 and 1. Thus it is important that Xd be at 
least of order unity, but larger values will yield little im¬ 
provement for the chosen detector geometry. We note 
Eq. [5] suggests Xd could be improved by several orders of 
magnitude by operation at the sub-Kelvin temperatures 
available in dilution refrigerators. However, as discussed 
in 0 , the susceptibility of practical paramagnetic insu¬ 
lators is well described by a Curie-Weiss relation of the 
form y = C/(T — Tew)-, where |Tc:vf| represents an effec¬ 
tive minimum temperature assuming the operating tem¬ 
perature is lower. Following Q, we use Tew = —2.1 K. 
We assume an operating temperature of 1.0 K, the lowest 
temperature at which the susceptibility obeys the Curie- 
Weiss law [4l|. This leads to an effective temperature of 
T = 3.1 K in Eq.j^and an estimate of Yrf = 0.53 (Table|T|, 
slightly below the value measured in ^ for single-crystal 
GGG. Optimization of other terms in Eq.[T0]is discussed 
below. 


TABLE I. Parameters used in the sensitivity computation. 


Parameter 

Value 

effective area of pickup coil, 

12.75 cm^ 

source mass density (BGO), pm 

7.13 g/cm^ 

detector spin density, Ud 

lO^Vcm® 

detector spin/Gd ion, s 

7/2 

detector susceptibility, Xd 

0.53 (SI) 

coupling efficiency, /3 

7.8 X 10"® 

sensor intrinsic noise, (pn. 

3/r4>o / V^Hz 

integration time, r 

10® s 


For more accurate estimates of the sensitivity of a prac¬ 
tical experiment to each potential in Eq. O we perform 
numerical calculations, in which all approximations used 
above are relaxed. The theoretical interactions from each 
potential due to a finite-sized source are computed at rep¬ 
resentative points in a hypothetical, practically-shaped 
detector and used to generate a spin excess profile as a 
function of the spatial coordinates. This profile is then 


TABLE 11. Dimensions of the detector and source masses used 
in the sensitivity calculations. 


Parameter 

Value 

detector width, Xd 

3.00 

cm 

detector length, yd 

3.00 

cm 

detector thickness, Zd 

0.76 

cm 

source mass width, Xm 

3.00 

cm 

source mass length, ym 

3.00 

cm 

source mass thickness, Zm 

1.00 

cm 

minimum gap, Zmin 

10 

fim 


used in a finite element (FEA) model to create a map 
of the induced flux in the region of the detector and 
pickup coil. 

The detector is broken up into subvolumes and the po¬ 
tentials calculated by Monte Carlo integration between 
the center point of each subvolume and the complete 
volume of the source mass. The detector and source 
dimensions (very similar to those used for the sample 
in 0 ) are given in Table IH] The induced spin orienta¬ 
tion at each subvolume location is obtained by repeating 
the calculation for many possible orientations and tak¬ 
ing the maximum. The integration models the modula¬ 
tion of the source mass (assumed sinusoidal) and records 
the results at several values of the separation (that is, 
the source phase) over a complete period. A minimum 
source-sample gap of 10 pm is used as in the analytical 
estimate. The maxima are then converted to magnetiza¬ 
tion vectors for each subvolume and phase using Eq. [T] 
and the parameters in Table m Fig. m shows a detector 
magnetization map for the case of the Vg+io interaction 
with A = 5 mm at a particular source phase. 

Magnetization maps are then entered into an electro¬ 
magnetic FEA model of the sample generated with the 
GOMSOL software package [i^. The model interpolates 
between the points of the input data to generate a com¬ 
plete magnetization profile within the detector volume 
and a complete field profile in the interior and exterior. 
An example field profile is shown in Fig. [2j 

From these prohles, the flux through the area of the 
pickup coils in the proposed experiment is calculated. 
Following 0 , we assume a planar gradiometer design for 
the pickup coil for the Vg+io and V12+13 interactions, to 
reduce common-mode backgrounds; the coil area in Ta¬ 
ble U matches that in 0 . For the V 4+5 interaction, the 
coil takes the planar “figure- 8 ” form shown in Fig. [2] for 
similar common-mode rejection, either sandwiched be¬ 
tween two halves of a split detector or wound through a 
small vertical hole in the center. For a particular interac¬ 
tion in Eq. [2] at a given range A, the experimental signal 
is calculated by taking the Fourier transform of the flux 
though the gradiometer coil as a function of the source 
phase, and scaling the result by the reduction factor in 
Eq. [S] For each value of A investigated, the source mass 
amplitude is optimized for maximum signal, resulting in 
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FIG. 2. (Color online) Top: map of detector mass magneti¬ 
zation induced by source in V 9+10 interaction. Vectors rep¬ 
resent magnetization (each an average of one of 343 subvol¬ 
umes). Bottom: magnetization (light arrows) and magnetic 
field (black arrows) induced by source in V4+5 interaction. 
Flux is calculated through cross section; dotted line shows 
placement of pickup coil for V4+5 experiment. In both cases 
A = 5 mm, and source (not shown) is slab of identical area 
centered above detector at instantaneous separation of 3 mm. 


amplitudes of order A. Finally, the SNR is obtained by 
dividing this result by the sensor intrinsic noise (Table|I|. 
The coupling constants in Eq.[ 5 ]are then adjusted so that 
SNR = 1 , resulting in the sensitivity curves in Fig. [31 

The velocity-dependent interactions V4+5 and V12+13 
are unconstrained for the case of polarized electrons. 
(We note that, if F4-1-5 is interpreted exclusively as an 
axial-vector interaction, the coupling i s tig htl y con¬ 
strained by electron spin-spin experiments [^, 1^.) 

For rough comparison, the bold solid line in the V4+5 plot 
is the limit on the corresponding coupling for polarized 
nucleons from the experiment at the Paul Scherrer Insti¬ 
tute ; the long-dashed line is the preliminary limit on 





FIG. 3. (Golor online) Projected sensitivity of proposed ex¬ 
periment to V 4 + 5 , V 9 + 10 , and Vi 2 -t-i 3 interactions, in which 
couplings are plotted versus range A (lower axes) or mass of 
unobserved boson (upper axes). Bold solid (dashed) curves 
are current (preliminary) limits, fine solid curves are indirect 
limits, dotted curves are proposals, dot-dash curve is Eq. ( 8 ] 
See text for explanations. 
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the nucleon coupling from [i^. Similarly, the bold solid 
line in the V 12+13 plot is the limit on the corresponding 
polarized nucleon coupling derived from the neutron spin 
rotation experiment at NIST 

The best limit on the 14+10 interaction for electrons 
derives from the axionlike particle torsion pendulum in 
the Eot-Wash group [ 2 ^. This limit is indicated by the 
bold line in the Vg +10 plot; the dotted line below it is the 
projected thermal limit. The sensitivity of the proposed 
experiment ranges from two to over ten orders of magni¬ 
tude greater in the range of interest. It is in good agree¬ 
ment with the analytical estimate from Eq. [5] at large 
A where A ~ suggesting that td better estimates the 
depth to which the detector is significantly magnetized. 

The fine solid lines in the V 4+5 and V 9+10 plots are the 
limits obtained from combining the constraints on gg or 
gp from the inferred cooling rate of SNI987A with the 
constraints on g^ from spin-independent short-range fifth 
force experiments . These apply to spin-0 interactions 
and are thus for illustrative purposes only in the case 
of Ej-i-s- The 14+10 limits are still more stringent than 
the projections for the direct search proposed, though 
the latter are more general and free from uncertainties 
related to dense nuclear matter effects in stars. 

The 14+10 plot also shows the limit inferred from the 
constraint on the electric dipole of the neutron {dn) [d^ . 
which sets bounds on gg via the QCD 0-term and thus is 
restricted to case of the (spin-0) axion. For generic light 
scalars unrelated to the strong CP problem, the limits 
from direct searches are more stringent than those in¬ 
ferred from EDM bounds over the range of interest [^ . 
Thus, correlating observations in EDM and macroscopic 
experiments could help distinguish axions from more 
generic light scalars. 

The light dotted lines in each plot are the projected 
sensitivity of the proposed experiment in [s^ . That pro¬ 
posal is sensitive to polarized electron coupling and com¬ 
petes directly, but our projections are stronger by about 
1-10 orders of magnitude for all interactions. A proposal 
for a direct experiment sensitive to the axion region in the 
14+10 plot is described in [5l|, though that experiment is 
sensitive to polarized nucleon coupling. 

Backgrounds .— The analysis above assumes the sen¬ 
sitivity of the experiment to be limited by the intrinsic 
noise of the SQUID sensor. Both GGG and BGO are 
good electrical insulators so Johnson noise levels are ex¬ 
pected to be very low, however, magnetic noise due to 
dissipation in the test masses is another possible statis¬ 
tical background. The spectral density of magnetization 
noise in the detector mass is given by [13, [111 : 




1 IkT^i" 
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where V is the volume and g" is the imaginary part of the 
complex permeability. While data on g" for the proposed 
GGG detector material are not available, it is expected to 
be too small (and even smaller in BGO, which has much 
lower susceptibility) for this noise to be observable. It 


was not observed, for example, in the experiment in [3, 
results from which were consistent with SQUID noise af¬ 
ter 5 days of integration time. 

Systematic effects are a greater concern. For example, 
the source can acquire a magnetization Mm = XmBo/g,o 
via the interaction of its susceptibility Xm with a stray 
field Bq. The changing flux through the pickup coil as 
the source oscillates above it can mimic a signal. BGO 
is diamagnetic with Xm = —1.9 x 10“^ at room temper¬ 
ature [53|, however there is evidence that it is weakly 
paramagnetic in low fields. Low-temperature data on 
Xm of pure BGO are not available, but the related com¬ 
pound a-Bi 203 exhibits Xm ~ 1-8 x I0“® at 4 K in low 
fields [13 . We estimate the flux as d(j) ~ {dB/dz)zoad, 
where dB/dz is the gradient above the center of a mag¬ 
netized disk of area Ud a distance Zd above the coil (Ta¬ 
ble HI, and zq is the largest source amplitude. Requiring 
d(j) < fpn/'/r leads to an upper limit on the allowed stray 
field of Bq Ri 10“^^ T. Assuming typical stray laboratory 
fields on the order of 10“^ T, the required shielding factor 
is 10®. This is quite modest compared to the experiment 
in [3 , which used two layers of superconducting lead foils 
and three layers of mu-metal sheets wound on frames to 
attain a measured shielding factor of 5 x 10^^ (which we 
assume here). 

Vibrations are potentially more troublesome. Motion 
of the detector-pickup coil assembly in phase with the 
source drive could induce a signal in the presence of a 
stray gradient dB^/dz, given by dcj) ~ {dBo/dz)Szdad, 
where Szd is the assembly vibration amplitude. This sig¬ 
nal is common-mode; the common-mode rejection ratio 
of the pickup coil used in Q was about 2 x 10^. Dividing 
the estimate for dcj) by this factor and again demanding 
d(f) < 4>nl\/T leads to a requirement on the stray gradient 
of dBojdz < 10“^^ T/m per micron of assembly vibra¬ 
tion. Assuming typical lab gradients of order I0“® T/m, 
this is within the shielding factor for vibrations less than 
about 100 p,m. If the vibrations cause the assembly to 
tilt in a stray field of 10“^^ T, the resulting flux signal 
estimate falls below the noise as long as the vibrations 
are less than about 30 fim. These effects can be studied 
by examining signals in the absence of a detector mass. 
The simple translating source mass assumed can be re¬ 
placed by a rotor with several segments of alternating 
density that pass over the detector at a multiple of the 
rotary drive frequency [5l|, [s^ , thus decoupling the drive 
frequency from the actual source modulation. 

Conclusions .— We have performed detailed calcula¬ 
tions of the projected sensitivity of an exotic interac¬ 
tion search with a paramagnetic insulating detector at 
cryogenic temperatures. The proposed technique affords 
the possibility to probe the interaction between macro¬ 
scopic test masses in near contact in a low-noise envi¬ 
ronment. Our results indicate either unique sensitivity 
to electron “monopole-dipole” interactions in the range 
below I mm, or improvements of more than ten orders 
of magnitude over existing experiments. The statistical 
limits in Fig. |3] represent the ultimate practical sensitiv- 
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ity of the experiment and are ambitious long-term goals. 
Results with reduced but competitive sensitivity, likely 
limited by systematics, are expected much sooner. The 
proposed technique is based largely on a proven design. 
Our primary purpose has been to show the potential sen¬ 
sitivity of that design, especially with the parameters of 
the existing detection scheme. Technical challenges asso¬ 
ciated with source mass translation in the cryostat and 
the related systematic backgrounds will certainly have 
to be addressed. Possible improvements include better 
SQUID coupling elficiency (/3 in Eq. [TUI), though this is a 


subtle optimization problem M- Increasing the detector 
area lEg.fTUlscales as y/od, given the near saturation of 
the demagnetization factor) is another possibility, though 
problems of test mass metrology and changes to SQUID 
coupling efficiency will warrant careful study. 

Acknowledgments .— The authors thank H. Gao, E. 
Smith, A. Holley, and W. M. Snow for useful discussions. 
This work is supported by National Science Foundation 
Grants PHY-1207656, PHY-1306942, Duke University, 
Los Alamos National Laboratory, and the Indiana Uni¬ 
versity Genter for Spacetime Symmetries (lUGSS). 


[1] E. Adelberger, J. Gundlach, B. Heckel, 

S. Hoedl, and S. Schlamminger, 

Prog. Part. Nucl. Phys. 62, 102 (2009) 

[2] I. Antoniadis, S. Baessler, M. Buch¬ 

ner, V. Fedorov, S. Hoedl, et al., 
Comptes Rendus Physique 12, 755 (2011) 

[3] J. Jaeckel and A. Ringwald, 

Ann. Rev. Nucl. Part. Sci. 60, 405 (2010) 

arXiv: 1002.0329 [hep-ph] 

[4] N. Arkani-Hamed, S. Dimopoulos, and 

G. Dvali, Phys. Lett. B429, 263 (1998) 

arXiv:hep-ph/9803315 [hep-ph] 

[5] S. Dimopoulos and A. A. Geraci, 
Phys. Rev. D68, 124021 (2003), 
arXiv:hep-ph/0306168 [hep-ph] 

[6] D. B. Kaplan and M. B. Wise, JHEP 0008, 037 (2000) 
arXiv:hep-ph/0008116 [hep-ph] 

[7] R. D. Peccei and H. R. Quinn, 
Phys. Rev. Lett. 38, 1440 (1977) 

[8] J. E. Moody and F. Wilczek, 
Phys. Rev. D30, 130 (1984) 

[9] Y. Kim, G.-Y. Liu, S. Lamoreaux, G. Visser, B. Kunkler, 
et al., arXiv:1104.4391 [nucl-ex] 

[10] A. O. Sushkov, S. Eckel, and S. K. Lamoreaux, 
Phys. Rev. A 79, 022118 (2009) 

[11] S. Eckel, A. O. Sushkov, and S. K. Lam¬ 
oreaux, Phys. Rev. Lett. 109, 193003 (2012) 

arXiv: 1208.4420 [physics.atom-ph] 

[12] G. Kittel and H. Kroemer, Thermal Physics, 2nd ed. (W. 

H. Freeman, 1980) p. 70. 

[13] A. Abragam, The Principles of Nuclear Magnetism (Ox¬ 
ford, 1978) p. 2. 

[14] F. L. Shapiro, Sov. Phys. Usp. 11, 345 (1968) 

[15] R. Bluhm and V. A. Kost- 

elecky, Phys. Rev. Lett. 84, 1381 (2000) 

arXiv:hep-ph/9912542 [hep-ph] 

[16] D. Budker, P. W. Graham, M. Ledbetter, S. Rajen- 
dran, and A. O. Sushkov, Phys.Rev. X4, 021030 (2014) 
arXiv:1306.6089 [hep-ph] 

[17] W.-T. Ni, S.-S. Pan, H.-C. Yeh, L.-S. Hon, and J.-L. 
Wan, Phys. Rev. Lett. 82, 2439 (1999) 

[18] A. N. Youdin, D. Krause, K. Jagannathan, L. R. Hunter, 
and S. K. Lamoreaux, Phys. Rev. Lett. 77, 2170 (1996) 

[19] M. Bulatowicz, R. Griffith, M. Larsen, J. Mirijanian, 

T. Walker, et al., Phys. Rev. Lett. Ill, 102001 (2013) 
arXiv:1301.5224 [physics.atom-ph] 

[20] G. Vasilakis, J. M. Brown, T. W. Kornack, and 
M. V. Romalis, Phys. Rev. Lett. 103, 261801 (2009) 


arXiv:0809.4700 [physics.atom-ph] 

[21] P. Chu, A. Dennis, G. Fu, H. Gao, R. Khati- 
wada, et al., Phys. Rev. D87, 011105 (2013) 
arXiv:1211.2644 [nucl-ex] 

[22] K. Tullney, F. Allmendinger, M. Burghoff, W. Heil, 
S. Karpuk, et al., Phys. Rev. Lett. Ill, 100801 (2013) 
arXiv: 1303.6612 [hep-ex] 

[23] R. C. Ritter, L. 1. Winkler, and G. T. Gillies, 

Phys. Rev. Lett. 70, 701 (1993) 

[24] G. D. Hammond, C. G. Speake, G. Trenkel, and A. P. 
Paton, Phys. Rev. Lett. 98, 081101 (2007) 

[25] B. R. Heckel, E. Adelberger, G. Gramer, T. Cook, 
S. Schlamminger, et al., Phys. Rev. D78, 092006 (2008) 
arXiv:0808.2673 [hep-ex] 

[26] S. A. Hoedl, F. Fleischer, E. G. Adelberger, and B. R. 
Heckel, Phys. Rev. Lett. 106, 041801 (2011) 

[27] S. Baessler, V. V. Nesvizhevsky, K. V. Protasov, 
and A. Y. Voronin, Phys. Rev. D75, 075006 (2007) 
arXiv:hep-ph/0610339 [hep-ph] 

[28] T. Jenke, G. Cronenberg, J. Burgdor- 
fer, L. Chizhova, P. Geltenbort, et al., 
Phys. Rev. Lett. 112, 151105 (2014) 

arXiv: 1404.4099 [gr-qc] 

[29] A. Serebrov, Phys. Lett. B680, 423 (2009) 

arXiv:0902.1056 [nucl-ex] 

[30] Y. Pokotilovski, Phys. Lett. B686, 114 (2010) 

arXiv:0902.1682 [nucl-ex] 

[31] A. K. Petukhov, G. Pignol, D. Jullien, and 

K. H. Andersen, Phys. Rev. Lett. 105, 170401 (2010) 
arXiv: 1009.3434 [physics.atom-ph] 

[32] C. B. Fu, T. R. Gentile, and W. M. Snow, 

Phys. Rev. D83, 031504 (2011) 

[33] A. K. Petukhov, G. Pignol, and 

R. Golub, Phys. Rev. D84, 058501 (2011) 

arXiv:1103.1770 [hep-ph] 

[34] B. A. Dobrescu and 1. Mocioiu, JHEP 0611, 005 (2006) 
arXiv:hep-ph/0605342 [hep-ph] 

[35] R. C. Ritter, C. E. Goldblum, W.-T. Ni, G. T. Gillies, 
and C. C. Speake, Phys. Rev. D 42, 977 (1990) 

[36] T. M. Leslie, E. Weisman, R. Khatiwada, and 
J. C. Long, Phys. Rev. D89, 114022 (2014) 
arXiv:1401.6730 [hep-ph] 

[37] S. K. Lamoreaux, Phys. Rev. A66, 022109 (2002) 
arXiv:nucl-ex/0109014 [nucl-ex] 

[38] C.-Y. Liu and S. Lamoreaux, 

Mod.Phys. Lett. A19, 1235 (2004) 

[39] S. Geller, in Physics of magnetic garnets, Vol. Proc. Inti. 
School Phys. “Enrico Fermi,” Course LXX, edited by 




S. Geller and A. Paoletti (Societta Italiana di Fisica, 
Bologna, 1978) p. 1. 

[40] M. Sato and Y. Ishii, 

Journal of Applied Physics 66, 983 (1989) 

[41] P. Schiffer, A. P. Ramirez, D. A. Huse, P. L. Gam- 
mel, U. Yaron, D. J. Bishop, and A. J. Valentino, 
Phys. Rev. Lett. 74, 2379 (1995) 

[42] Comsol, Version 4-4 i Comsol, Inc., 1 New England Ex¬ 
ecutive Park, Burlington, MA 01803 USA. 

[43] B. R. Heckel, W. A. Terrano, and E. G. Adelberger, 
Phys. Rev. Lett. Ill, 151802 (2013) 

[44] S. Kotler, R. Ozeri, and D. F. J. Kimball, 
arXiv:1501.07891 [physics.atom-ph] 

[45] F. M. Piegsa and G. Pignol, 

Phys. Rev. Lett. 108, 181801 (2012) 
arXiv:1205.0340 [hep-ex] 

[46] Y. Zhang, G. Sun, S. Peng, G. Fu, H. Guo, et al, 
arXiv:1412.8155 [nucl-ex] 

[47] H. Yan and W. M. Snow, 

Phys. Rev. Lett. 110, 082003 (2013) 


arXiv:1211.6523 [nucl-ex] 

[48] G. Raffelt, Phys. Rev. D86, 015001 (2012) 

arXiv: 1205.1776 [hep-ph] 

[49] G. Baker, D. Doyle, P. Geltenbort, K. Green, M. van der 
Grinten, et al., Phys. Rev. Lett. 97, 131801 (2006) 
arXiv:hep-ex/0602020 [hep-ex] 

[50] S. Mantry, M. Pitschmann, and M. J. 

Ramsey-Musolf, Phys. Rev. D90, 054016 (2014) 

arXiv:1401.7339 [hep-ph] 

[51] A. Arvanitaki and A. A. Geraci, 
Phys. Rev. Lett. 113, 161801 (2014) 

arXiv: 1403.1290 [hep-ph] 

[52] S. Eckel, A. O. Sushkov, and S. K. Lamoreaux, 
Phys. Rev. B 79, 014422 (2009) 

[53] S. Yamamoto, K. Kuroda, and M. Senda, 
IEEE Trans. Nucl. Sci. 50, 1683 (2003) 

[54] E. Kravchenko, V. Orlov, and M. Shlykov, Russ. Chem. 
Rev. 75, 77 (2006). 

[55] Y. J. Chen, W. Tham, D. Krause, D. Lopez, E. Fis- 
chbach, et al., arXiv:1410.7267 [hep-ex] 


